Reflection of light and heavy holes from a linear potential barrier. 
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In this paper we study reflection of holes in direct-band semiconductors from a linear potential 
barrier. It is shown that light-heavy hole transformation matrix depends only on dimensionless 
product of the light hole longitudinal momentum and the characteristic length determined by the 
slope of the potential and it doesn't depend on the ratio of light and heavy hole masses, provided this 
ratio is small. This coefficient is shown to vanish both in the limit of small and large longitudinal 
momenta, however the phase of a reflected hole is different in these limits. An approximate analytical 
expression for the light-heavy hole transformation coefficient is found. 

PACS numbers: 7155.Eq, 7210.Fk 



I. INTRODUCTION 



Direct band A3B5 semiconductors are very challenging 
both for theoretical and experimental applications. In 
particular, heterostructures based on these compounds 
are widely studied. There have been published many pa- 
pers devoted to the boundary condition problem in het- 
erostructures J|, ||, U, [| @) 01- These articles mainly 
dealt with abrupt interface potentials, where appropri- 
ate boundary conditions were sought to relate the wave 
functions at the different sides of the interface. The most 
common approach used to find the boundary conditions 
is to integrate envelope function equations through the 
interface [g, ||, ||, [| . In this way, it appears that there is 
no mixing between light and heavy holes at normal inci- 
dence [[|. However, there is no general reason justifying 
the integration procedure. Moreover, the tight binding 
model and some numerical calculations (see M and refer- 
ences therein) are in direct contradiction with the above 
statement. In particular, for a GaAs/AlAs interface the 
light-heavy hole transformation coefficient for the nor- 
mal incident flux was found to be of order unity || . In a 
purely phenomenological analysis, relying on group theo- 
retical formalism Q , the boundary conditions were found 
to be determined by a number of phenomenological pa- 
rameters which should be found either from microscopical 
calculations or from fitting to experimental data. How- 
ever, this approach hardly allows following the dynamics 
of the reflection matrix as a function of energy and lateral 
momentum of incident holes. On the other hand there 
is possible a completely different situation wherein the 
potential near the turning point is smooth. The bound- 
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ary of this type may appear either in a compound with 
relatively slow variation of composition or in a structure 
subjected to an external electric field. Smoothness of the 
potential justifies use of the envelope function equations 
and hence eliminates boundary condition uncertainties. 

The formulation of the problem suggests that the semi- 
classical approach is to be used. However, being applied 
to Luttinger equations, it gives new essential features as 
compared to standard WKB formulas. The situation of 
a smooth linear potential was first considered in for 
a gapless system, and later, in fLo| , for the case of A 3 B 5 
semiconductors. In the latter paper, it was shown that 
there is a light-heavy hole transformation even in a slow 
varying potential if the incident particle flux is not nor- 
mal to the interface. In their calculation, the authors 
used an infinitely large heavy hole mass, which resulted 
in divergence of the wavefunction in a certain underbar- 
rier point at which the kinetic energy of a hole was equal 
to zero. In the present paper we take the ratio of heavy 
and light hole masses to be large but finite and avoid di- 
vergences of this type. We note finally that though we 
assume that the hole masses and spin-orbital splitting 
don't depend on the coordinate, their slow linear varia- 
tions can be easily included into the model. 



II. SOLUTION OF LUTTINGER EQUATIONS 
IN A LINEAR EXTERNAL POTENTIAL 



Let us start from the conventional 8x8 Hamiltonian 
in a spherical band approximation, i.e. assuming that 
the Luttinger parameters 72 and 73 are equal to each 
other This assumption is usually well justified 

in the vicinity of T-point. Instead of total angular mo- 
mentum basis, we use another one || |l2| 

| s t>, \si), M>, M>, M>, \vi), M>, H>, 

(i) 

which allows us to write the equations for electron and 
hole wave functions in a compact analytical form: 
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(E g + <j)(x) - E)^ s - 67V* = 

(-E + <P(x) - ^) * - £ 7 W S + -^(f! + 4f 2 )V(W) - -^(71 - 27~ 2 ) V x [V x *] + i^[a X *] = 0, (2) 
\ o J 2m 2m 6 

I 



where 71 and 72 are generalized Luttinger parame- 
ters |n| , 7 is the Kane parameter having dimension of ve- 
locity, m is a free electron mass, a = (a x , <J y ,a z ) are Pauli 
matrices, E g and A so are the band gap and spin-orbital 
splitting, <p(x) is the external potential which will be con- 
sidered hereafter to be a linear function of x-coordinatc: 
4>(x) = ax. We use the units where Plank constant is 
equal to unity. The equation above is a generalization 
of that suggested in (^] to the case of a finite heavy hole 
mass. Obviously the wave function is to be found as a 



plain wave propagating along the lateral directions mul- 
tiplied by a function of x coordinate: 

$r A (r) = tp x (x) jiyy+i*** 

with A = s, x, y, z. Note that ^ \ and ip\ are the spinors. 
Without loss of generality choose q z = 0. We further 
simplify our calculations assuming that E g and A so are 
much larger than E and cj>. This corresponds to the 4x4 
model 0). Then (|^) becomes: 



mih +3m h d _JU +2aX -2E 

Am h m lh dx 2 m h 



±i 



J2_^_ ± 3c(m h -m lh ) d + _ E 
2m h Am h m lh dx 



2m h 



3q(m h —mi h ) d 
4m h mi h dx 



ax — E 



1 d 2 _ q (rni h +3m h ) 
m h dx 2 im h m lh 



2ax - 2E 



\ ( \ 

1px(x) 



= 



(3) 



where m/, and mih are heavy and light hole effective 
masses related to Luttinger parameters as: 

m m 4to7 2 

=71-272 = — — +71 +2 72 . 

m h mih ZE g 

The upper sign in (g) corresponds to the spinor satisfy- 
ing: 

<Jz^x = Ipx, °z^y = Ipy, zi> z = ~1pz, 

while the lower does to: 

cr z i> x = -ip x , <J z ip v = -tpy, a z ip z = t/jg, 

Obviously the two sets of solutions are obtained form 
each other via the time reversal transformation: q — > —q 
and ip —* ip*. Therefore we will consider only the system 
(0) with the upper sign. 

It proves to be convenient to introduce new linear com- 
binations: 

ip + = ip x + 2iip y , ip- = -2ipx - i?Py (4) 
Then the equation (J3J) can be written in an elegant way: 

5 - q 2 ^ + ) - 2\iX^+ =qF 
5 (ipl - q 2 i>-) - 2A ; a;V- = F' , (5) 

ip' + + 2t/j'_ - q{2ip+ +ip-) = -3F(1 - S)- 1 



where 



mi 



« 1, 
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F is an unknown function to be determined and intro- 
duced for convenience. After rescaling coordinates and 
momenta: 



A, 



-1/3 



A, l /» 



x = x- 



q = r] 



1-5 



F = F- A ' 



1/3 



1 - 5 



1-5 

equation (||) becomes: 

5 O" - rfip+) - 2(1 - 6)xip+ = T]F 
5 O" - rfip-) - 2(1 - 5)xip- = F' 
tp' + + 2tp'_ - r/(2^+ + ip-) = -3F(1 - 5)- 1 



(6) 



It is convenient to perform Fourier transform and then 
substitute: 



= ijj+(k) exp (i^jz 



5) 



ip-(k) = ip-(k)exp ^2(T=3) 
F(k) = G{k) exp 



(7) 



2(1-5) 



Finally, we obtain a second order differential equation for 
the function G: 



d 2 G 
dk 2 



.k 2 + r] 2 dG 
' 3 ~dk 



G = 0. 



(8) 
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The boundary condition corresponding to the incident 
light hole is defined as the asymptotic at k — > oo : 



The other asymptotic at k 
G -> Ar H vV + k 2 



i exp 



A.- 3 



-oo is: 



< 3 tan » exp ( ~ 



rfk 



1.3 



(9) 



{k 2 +T] 2 f /2 '' 



where ru and r^; are the light-light and light-heavy hole 
reflection coefficients. Below we relate them to the light 
and heavy hole fluxes in real space. Similar expressions 
can be written for the heavy holes to define Thh and rih'- 



G^B 



{k 2 +T] 2 ) 3/2 



at k 



oo 



(11) 



and 

G- 



Brih^Jrf + k 2 



-Br hh 



tan- 1 i I % 

( 2 tan i exp ( - 



(k 2 + v 2 ) 3/2 



at k — ► — c 



(12) 



Time reversal symmetry of (jg) gives simple relations: 



\m\ +rf ll n h 

r u r hi + r* hl r hh = 



or rr* = I, 



(13) 



where / is a 2 x 2 unit matrix. Note that r is not unitary 
and it depends only on a single parameter r\. 

Let us proceed with the analysis of (^) for the case of 
large r\. It appears that in this limit solving this equation 
is similar to finding a superbarrier reflection coefficient in 
ordinary quantum mechanics |l3| . This analogy suggest 
shifting the variable k from the real axis. For positive r\ 
we substitute: 



yielding: 



d 2 G X 



k = ir) + £ 



dG 



5rj 



W -^ + ^-^-^)G = 0. (14) 

The boundary condition corresponding to the incident 
light hole is: 

G^^ 1/4 (2^ + 3/4 e"^" 3+i ^"^ a atf^oo. (15) 
Provided ( < ?j, ( |i"4| ) simplifies to: 
d 2 G 2r]£ dG 5ry 



de ' ~^ + y g -°- 



(16) 



The solution to this equation, which is formally can be 
expressed through Whittaker functions, is most conve- 
niently given via the integral representation: 



G = Ae-^\ 



33/4 
\/2tt 



oo+"i0 



J t L '*e-%e it *dt. (17) 



-00-f iO 



The particular choice of a multiplier before the inte- 
gral ensures coincidence of ( |l7| ) and (|l~J) in the region 
yj rj -1 -C £ -C rj. This can be verified by applying 



(10) the steepest descend integration to (17). In the region 



-77 < £ < 



'rj 1 we find asymptotics corresponding 



to rhi and ru. However, as noted above the ru coefficient 
must be set to be 1 and not the value obtained in our ma- 
nipulation. In fact loss of accuracy for some branches of a 
solution is quite common when asymptotical expressions 
arc extended into the complex plain. On the other hand 
the rhi coefficient is found with exponential accuracy: 

r w « e-i"V*T ^(2r?) 7 / 4 r(l/4), (18) 

where T denotes the T-function. If 77 is large and negative, 
one has to shift variable to k — — + Coming through 
the same steps we can get: 



r h i « e gV e i 4 



3 9/4 



2 7 /40f 



|»?r/ 4 r(3/4). (19) 



In a similar fashion we can provide analysis for small 
values of 77. As will be clear later it is sufficient to restrict 
calculations only to the case 77 = 0. Then (|J) has two 
simple solutions: 



G 2 (k) = 3ke l - +l ^P 



t 2 



(20) 
(21) 



where P stands for the principal value of the integral. 
The coefficients before G\(k) and G2(k) are chosen in 
the way that at k — + 00 these functions coincide with the 
defined above asymptotics of the light and heavy holes 
respectively. It is not difficult to observe that 



ru = i 



ru = 



rih = -^r(2/3) r hh 



(22) 



This matrix obviously satisfies conditions (P~q). 

Now let us return to real space. From (||) and (Q) we 
obtain 

00 00 
ip + (x) = 77 J dfc e 4 5<i^)[ i #+^fc]+^y* G(t)dt, (23) 

— OO k 
OO OO 

i/i-(x) = i [ dke^^^ Wkl \ +lki [ tG(t)dt. (24) 



Integrating equation (Jg) itself and multiplied by k, one The heavy hole component is given by 
hnds: 



' Gm= 2b -«fe) f g(fc)e -*(*^»)V+i* 



k 

OO 



tG {t)dt = 2(2 ^ + fc) [G(k)e : 
rj 



i^G(k){25) 
2 



+^ fe ) \ _^G(fc). (2$(*) 
* »7 



In the limit \rj\ » 1, rjj is close to 1 and using (|J), (25), 
and (^6|) we can find that the light hole component of the 
wave function is approximately equal to: 



V°(x) 
^_(x) 
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2^ 
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t28) 



iAr h i 



3v37r?7 



rje 2 



- OO 





3v37r?7 
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; „A tan -1 - 
fee 2 *> 

Vfc 2 + ?7 2 



1 '2.-16 3 



' 2-2rt 3 



'(29) 



(30) 



We used x — y basis for the last equations since it is 
more suitable for heavy holes. From (^) and (^8|) we 
deduce semiclassical expressions for the light hole wave- 
functions valid at large negative x, provided \r)\ 1. 



3 f tA , 71 ./-i^-cos f fc,(x')cfi' + i tan- 1 Mf) - 



. . sin f fe ; (x')c/i' + i tan- 1 M*l 



Here a; is a classical turning point and fc;(x) is a position-dependent light hole wavevector: 

2 1 - S/2 



r 



(31) 



Similarly we can deduce semiclassical wave functions 
valid for i/ < 1. The only difference, which appears as 
compared to ( |3l| ) is the absence of (tan- 1 /?;/?]) in the 
argument of cosine. 

We note that at 5 — > expressions (^9|) and (^0|) di- 
verge at x — > 0. In fact, as was noted in jn|, this phe- 
nomenon is closely related to the divergence of the elec- 
tric field in electromagnetic waves incident to the media 
with linearly increasing dielectric constant |l4[ . It can 
be shown that the system (||) is equivalent to the fourth 
order differential equation with a small coefficient at a 
higher derivative. And the equations under considera- 
tion have precisely the form of the Maxwell equations for 
the electric and magnetic fields. The parameter 5 is anal- 



ogous to the coupling between plasmons and the electro- 
magnetic field. In other words, the effect of heavy holes 
on the light hole wave functions is similar to account 
of spatial dependence of dielectric constant in plasma. 
However, in order to get a sharp resonance, 8 should be 
very small. For real semiconductors its value is usually 
bounded between 10 -2 and 10 , which appears to be 
insufficient to observe the peak. On the contrary, for 
plasma the corresponding quantity (S p ~ T/mc 2 ) is sev- 
eral orders of magnitude less [Q, and hence the reso- 
nance can exist. 

In a similar manner WKB wave functions for the heavy 
holes can be found. Thus for \rj\ > 1 in a classically 
allowed region we get: 



5 



where 



2,^/2 2(1-5). 

Gh = ~ 2(l-S) r] ' = " ~ V * X 

r 



Note that the turning point for heavy holes (an) is closer 
to the origin then that for the light holes as should be 
expected. 

The next step in our analysis is to relate the r-matrix 
in momentum space defined in (|lo| ) to the _R-matrix in 



real space. To do this we have to define the basis func- 
tions for the incident and reflected waves. Like we have 
observed already, it is convenient to work in the +, — 
representation for the light holes and x, y representation 
for the heavy holes. 



Pi 



1 



=pz f ki(x')dx'^f % tan 



-l fcj (g) i_ 



kh(x) 



c 



Ti f k h (£')d£'±% tan" 1 ^^±if 



(33) 



(34) 



where <^ refer to the incoming and outgoing waves. The x-components of the fluxes associated with these basis 
wavefunctions are equal to: 



mi 
h 

rrih 



i dip* , 



dx 

d^lx 
dx 



±- 



mi 



qvh,xv h ,y + 2 ^ h -y~dx~ 



±- 



(35) 
(36) 



Here we explicitly inserted A/ and ?L The choice of the 
basis functions ( |33] ) and (|4|) is now evident, since they 
carry out equal fluxes. Now it is not hard to relate the 
R and r matricies: 
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Rlh 




i Rhl 


Rhh j 
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ru 



_3\/3rj_ 



2y/l-S/2 



rid 



2y/l-S/2 

rhh 



rih 



(37) 



The unitarity condition for R requires that 
, ., . 4 



\ni\ 



27^ 



(38) 



where we neglected by <5 compared to 1. All other rela- 
tions are just a consequence of (fL3h, in particular: 



Rhi\ = \R, 



lh 



\Rhh \ — \Rll\- 



Using ([L8|) , ( jl9| ) and fl37| ) we can find the approximate 
expression for the light-heavy hole transformation coeffi- 



cient | Rm | 2 valid for large r\ 
,2 1 2 %/2 



|iZ«|' 



|i?w| 2 



if /2 T 2 (l/i)e~^ at 77 -> 00, (39) 



13^3 1 ^ 



1^(3/4) at ^ -00, (40) 



7T 2\/2 V 3/2 
At small 77, from (E2) and (B7|) we obtain: 



\Rhi\ 2 = \Rih\ 2 ~ ^^r 2 (2/3) at 77^0. (41) 



Using asymptotics (|39j) , (|40j) and (|4lj) as well as the exact 
numerical results it is possible to find an interpolating 
function which gives an excellent description of |-R/j.;| 2 
for all values of ij (see Figure 1): 



^ ^1^(2/3) e 1 - 3 ^ 



,- 9 31/^(2/3) , a,, . 1 
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at 77 > 0, (42) 
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FIG. 1: Dependence of the light-heavy hole transformation 
coefficient \Rhi\ 2 on 77, obtained from numerical solution of 
(^) - solid line, and the approximate expressions (^2|), (^) - 
dotted line. The insert shows the relative phase of the light 
hole reflection coefficient Ru. 



FIG. 2: Real and imaginary part of %))+ as a function of x- 
coordinate for the case of the incident light hole. The top 
graphs refer to the large mixing between the hole branches 
(77 = 1) and the bottom graphs do to the small mixing (77 = 
2.5). 



tial. In particular one finds: 



III. DISCUSSION 



From the expression (|43|) it follows that the light-heavy 
hole transformation coefficient vanishes both in the limit 
of small and large 77. The maximums occur at 77 close 
to ±1, see Figure 1. Note that the maximum reflection 
coefficient at positive 77 is much larger than that at neg- 
ative. An important result is that the dependence of the 
reflection matrix R (and transformation coefficient Rm 
in particular) on the dimensionless longitudinal momen- 
tum 77 is not sensitive to the ratio of the light and heavy 
hole masses. Also the potential slope completely scales 
out and the crossover from small to large q occurs even 
in the case of a very slow varying potential. 

In the insert of Fig. 1 we plot the phase of the light hole 
reflection coefficient. It is equal to the conventional value 
for large 77. The discontinuity of the phase at 77 = is 
due to inclusion of ± tarT 1 (k/rj) into the basis functions 
([33|) and (|34|). Similar curve but with the opposite sign 
is valid for the heavy holes. 

Now let us compare transformation coefficient studied 
in this paper with that obtained for a scattering on an 
abrupt potential. There is no unique and justified ap- 
proach to derive the boundary conditions for envelope 
functions However, for the infinite barrier the most 
natural and used requirement is the vanishing of wave 
function components at the interface (see for example 
refs ||, |ll|). Using these boundary conditions and 
basis functions ( p3| ) and ( |34| ) with ki and kh being con- 
stants independent of a;, it is not hard to derive reflection 
matrix for the scattering from an infinite abrupt poten- 



\K h \ 2 



\r; 



a |2 
lh\ 



I2kik h q 2 



a |2 



4(k l k h + q 2 ) 2 + q 2 (k h -k l ) 2 
Aihkh-q 2 ) 2 +q 2 (k h + k l ) 2 



4(kik h + q 2 ) 2 + q 2 {k h - h) 2 
Introducing the incidence angle for the light holes: 
• ± 1 



,(44) 



(45) 



kf + q 2 



and using i C 1 we get: 



\R a hl \ 2 = 



6v6<5sin (Acos< 



(sin* 



4 cos 2 



0(5). 



(46) 



In a more general case of nonspherical band (72 7^ 73) 
light and heavy hole transformation coefficients were nu- 
merically investigated in ||. Expression (46), contrary 
to ([l3]), has explicit dependence on 6. It vanishes in both 
limits of normal and lateral incident light hole. How- 
ever, as we mentioned above, use of vanishing at inter- 
face boundary conditions is controversial and needs fur- 
ther justification. Note that the reflection matrix from a 
linear smooth potential is found without making any spe- 
cific assumptions about semiconductors. So the result of 
the present analysis is quite general and doesn't depend 
on microscopical details. 

Figure 2 shows the wave functions for the incident light 
hole for two different values of 77 corresponding to large 
and small mixing between the two hole branches. High 
frequency oscillations correspond to the heavy hole com- 
ponent, which is obviously essential for 77 = 1 and rel- 
atively small for 77 = 2.5. Note that as was mentioned 
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earlier at 5 — > there appears a sharp resonance in the 
wave function near x — 0. However it is hardly seen in 
the graphs. The reason is that the height of the peak, as 
it follows from (f|, |o|), is of the order (5" 1 / 3 exp(-2/9 T? 3 ) 
and the peak itself exists only when 3> 1. Real values 
of S for the semiconductors are not less then 10~ 2 , which 
is insufficient for observing the resonance. 

In conclusion, the light-heavy hole transformation co- 
efficient for reflection from a linear potential is found to 
be a function of a dimensionless parameter 77, which is 
proportional to the longitudinal momentum of the light 
holes, and to be independent of the ratio of the light and 
heavy hole masses if this ratio is small. This function 
vanishes in both limits of large and small rj, however the 



phases of light-light and heavy-heavy hole reflection co- 
efficients are different in these limits. An approximate 
analytical expression for the transformation coefficient is 
found (see (fl3|)). Account of the finite heavy hole mass 
is shown to be responsible for disappearance of an un- 
derbarrier resonance in the light hole wave function pre- 
dicted in (hJ. 

After the paper was written authors became aware of 
the work [ ft6| , where some of the results, in particular 
asymptotical expressions ( |39| ) and ( fio| ) , were already ob- 
tained. 

This work was partially supported by the Russian 
Foundation for Basic Research. One of the authors is 
grateful to I. Yu. Solov'cv for useful discussions. 
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